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SUPERCONDUCTING QUBITS AND THE PHYSICS OF 
JOSEPHSON JUNCTIONS 



John M. Martinis , Kevin Osborne 



1 Introduction 

Josephson junctions are good candidates for the construction of quantum 
bits (qubits) for a quantum computer [1] . This system is attractive because 
the low dissipation inherent to superconductors make possible, in principle, 
long coherence times. In addition, because complex superconducting cir- 
cuits can be microfabricated using intcgrated-circuit processing techniques, 
scaling to a large number of qubits should be relatively straightforward. 
Given the initial success of several types of Josephson qubits [2-10], a ques- 
tion naturally arises: what are the essential components that must be tested, 
understood, and improved for eventual construction of a Josephson quan- 
tum computer? 

In this paper we focus on the physics of the Josephson junction because, 
being nonlinear, it is the fundamental circuit element that is needed for 
the appearance of usable qubit states. In contrast, linear circuit elements 
such as capacitors and inductors can form low-dissipation superconducting 
resonators, but arc unusable for qubits because the energy-level spacings 
are degenerate. The nonlincarity of the Josephson inductance breaks the 
degeneracy of the energy level spacings, allowing dynamics of the system 
to be restricted to only the two qubit states. The Josephson junction is 
a remarkable nonlinear element because it combines negligible dissipation 
with extremely large nonlincarity - the change of the qubit state by only 
one photon in energy can modify the junction inductance by order unity! 

Most theoretical and experimental investigations with Josephson qubits 
assume perfect junction behavior. Is such an assumption valid? Recent 
experiments by our group indicate that coherence is limited by microwave- 
frequency fluctuations in the critical current of the junction [10]. A deeper 
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understanding of the junction physics is thus needed so that nonideal behav- 
ior can be more readily identified, understood, and eliminated. Although 
we will not discuss specific imperfections of junctions in this paper, we want 
to describe a clear and precise model of the Josephson junction that can 
give an intuitive understanding of the Josephson effect. This is especially 
needed since textbooks do not typically derive the Josephson effect from 
a microscopic viewpoint. As standard calculations use only perturbation 
theory, we will also need to introduce an exact description of the Josephson 
effect via the mesoscopic theory of quasiparticle bound-states. 

The outline of the paper is as follows. We first describe in Sec. 0thc 
nonlinear Josephson inductance. In Sec. |3|we discuss the three types of 
qubit circuits, and show how these circuits use this nonlinearity in unique 
manners. We then give a brief derivation of the BCS theory in Sec. 0J 
highlighting the appearance of the macroscopic phase parameter. The 
Josephson equations are derived in Sec. El using standard first and second 
order perturbation theory that describe quasiparticle and Cooper-pair tun- 
neling. An exact calculation of the Josephson effect then follows in Sec. |S] 
using the quasiparticle bound-state theory. Section [71 expands upon this 
theory and describes quasiparticle excitations as transitions from the ground 
to excited bound states from nonadiabatic changes in the bias. Although 
quasiparticle current is typically calculated only for a constant DC voltage, 
the advantage to this approach is seen in Sec. HO where we qualitatively 
describe quasiparticle tunneling with AC voltage excitations, as appropri- 
ate for the qubit state. This section describes how the Josephson qubit 
is typically insensitive to quasiparticle damping, even to the extent that a 
phase qubit can be constructed from microbridge junctions. 

2 The Nonlinear Josephson Inductance 



A Josephson tunnel junction is formed by separating two superconducting 
electrodes with an insulator thin enough so that electrons can quantum- 
mechanically tunnel through the barrier, as illustrated in Fig. ^ . The 
Josephson effect describes the supercurrent Ij that flows through the junc- 
tion according to the classical equations 



where $o = h/2e is the superconducting flux quantum, Iq is the critical- 
current parameter of the junction, and 5 = <p L ~ (f> R and V are respectively 
the superconducting phase difference and voltage across the junction. The 



/,/ = I sin S 

v = i^— 

2?r dt ' 



(2.1a) 



(2.1b) 
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Fig. 1. Schematic diagram of a Josephson junction connected to a bias voltage 
V. The Josephson current is given by Ij = Josin<5, where S — cj> L — <j> R is the 
difterence in the superconducting phase across the junction. 



dynamical behavior of these two equations can be understood by first dif- 
ferentiating Ea. l2.lal and replacing d5 / dt with V according to Ea. l2.lbl 

dlj , 2tt , , 

- = I co S 6-V. (2.2) 

With dlj/dt proportional to V, this equation describes an inductor. By 
defining a Josephson inductance Lj according to the conventional definition 
V = Ljdlj/dt, one finds 

Lj = - ^ . (2.3a) 
2tt1o cos o 

The 1/ cos S term reveals that this inductance is nonlinear. It becomes 
large as 6 — > tt/2, and is negative for ir/2 < S < 3tt/2. The inductance at 
zero bias is Ljo — <i>o/27rio- 

An inductance describes an energy-conserving circuit element. The 
energy stored in the junction is given by 

Uj = J IjVdt (2.4a) 
7 sin<5/-> (2.4b) 

27T dt 



2ir 



sinSdd (2.4c) 
cos 6. (2.4d) 



This calculation of energy can be generalized for other nondissipative circuit 
elements. For example, a similar calculation for a current bias gives L^ias = 
— (I$o/2tt)5. Conversely, if a circuit element has an energy U(S), then the 
current-phase relationship of the element, analogous to Ea. l2.lal is 

27TdU(S) 

Ij{6) = ■ (2 - 5) 
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A generalized Josephson inductance can be also be found from the second 
derivative of U , 

1 /2vr\ 2 d 2 U{6) 



Lj VW dS 



r2 



(2.6) 



The classical and quantum behavior of a particular circuit is described by 
a Hamiltonian, which of course depends on the exact circuit configuration. 
The procedure for writing down a Hamiltonian for an arbitrary circuit has 
been described in detail in a prior publication [11]. The general form of the 
Hamiltonian for the Josephson effect is Hj = Uj. 

3 Phase, Flux, and Charge Qubits 



A Josephson qubit can be understood as a nonlinear resonator formed from 
the Josephson inductance and its junction capacitance. nonlinearity is 
crucial because the system has many energy levels, but the operating space 
of the qubit must be restricted to only the two lowest states. The system 
is effectively a two-state system [12] only if the frequency uj\q that drives 
transitions between the qubit states < — > 1 is different from the frequency 
u)2i for transitions 1 < — > 2. 

We review here three different ways that these nonlinear resonators can 
be made, and which are named as phase, flux, or charge qubits. 

The circuit for the phase-qubit circuit is drawn in Fig. Ufa). Its Hamil- 
tonian is 

1 ^ 2 I <P 7 I$ ~ ,„ ,s 

H= 2C Q (3 - 1} 
where C is the capacitance of the tunnel junction. A similar circuit is 
drawn for the flux-qubit circuit in Fig. Efb), and its Hamiltonian is 

The charge qubit has a Hamiltonian similar to that in Eq. 13.11 and is 
described elsewhere in this publication. Here we have explicitly used no- 
tation appropriate for a quantum description, with operators charge Q and 
phase difference S that obey a commutation relationship [5, Q] — lei. Note 
that the phase and flux qubit Hamiltonians are equivalent for L — > oo and 
I = <f>/L, which corresponds to a current bias created from an inductor with 
infinite impedance. 

The commutation relationship between 8 and Q imply that these quan- 
tities must be described by a wavefunction. The characteristic widths of 
this wavefunction are controlled by the energy scales of the system, the 
charging energy of the junction Ec — e 2 /2C and the Josephson energy 
Ej = 7o$o/27r. When the energy of the junction dominates, Ej ^> Ec, 
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Fig. 2. Comparison of the phase (a), flux (b), and charge (c) qubits. Top row 
illustrates the circuits, with each "X" symbol representing a Josephson juncton. 
Middle row has a plot of the Hamiltonian potential (thick line), showing qualita- 
tively different shapes for three qubit types. Ground-state wavefunction is also 
indicated (thin line) . Key circuit parameters are listed in next row. Lowest row 
indicates variations on the basic circuit, as discussed in text. The lowest three 
energy levels are illustrated for the phase qubit (dotted lines). 

then 5 can almost be described classically and the width of its wavefunction 
is small (5 — (5) ) -C 1- In contrast, the uncertainty in charge is large 
(Q 2 - (Q) 2 ) » (2e) 2 . 

If the Josephson inductance is constant over the width of the 5 wave- 
function, then a circuit is well described as a Lj-C harmonic oscillator, and 
the qubit states are degenerate and not usable. Usable states are created 
only when the Josephson inductance changes over the S- wavefunction. 

The most straightforward way for the wavefunction to be affected by the 
Josephson nonlinearity is for 5 to have a large width , which occurs when 
Ej ~ Eq- A practical implementation of this circuit is illustrated in Fig 
Etc) , where a double-junction Coulomb blockade device is used instead of 
a single junction to isolate dissipation from the leads [2,4]. Because the 
wavefunction extends over most of the -cos 5 Hamiltonian, the transition 
frequency uj\o can differ from uj2\ by more than 10 %, creating usable qubit 
states [13]. 

Josephson qubits are possible even when Ej ^> Ec, provided that the 
junction is biased to take advantage of its strong nonlinearity. A good 
example is the phase qubit [6], where typically Ej ~ 10 4 Ec, but which 
is biased near 8 ;$ tt/2 so that the inductance changes rapidly with S (see 
Ea. l2.3a(l . Under these conditions the potential can be accurately described 
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by a cubic potential, with the barrier height AU — > as I — ► Iq. Typically 
the bias current is adjusted so that the number of energy levels in the well is 
~ 3 — 5, which causes u>iq to differ from uj 2 i by an acceptably large amount 
~ 5 %. 

Implementing the phase qubit is challenging because a current bias is 
required with large impedance. This impedance requirement can be met by 
biasing the junction with flux through a superconducting loop with a large 
loop inductance L, as discussed previously and drawn in Fig. Ufa). To 
form multiple stable flux states and a cubic potential, the loop inductance 
L must be chosen such that L > 2Lj Q . We have found that a design 
with L ~ 4.5Lj is a good choice since the potential well then contains the 
desired cubic potential and only one flux state into which the system can 
tunnel, simplifying operation. 

The flux qubit is designed with L j§ Ljq and biased in flux so that 
(S) = 7r. Under these conditions the Josephson inductance is negative and 
is almost canceled out by L. The small net negative inductance near 5 = n 
turns positive away from this value because of the 1/ cos S nonlinearity, 
so that the final potential shape is quartic, as shown in Fig. ^b). An 
advantage of the flux qubit is a large net nonlinearity, so that wio can differ 
from W21 by over 100 %. 

The need to closely tune L with Ljq has inspired the invention of several 
variations to the simple flux-qubit circuit, as illustrated in Fig. Hfb). One 
method is to use small area junctions [7] with Ej ~ lO-Ec, producing a large 
width in the 5 wavefunction and relaxing the requirement of close tuning of 
L with Ljq. Another method is to make the qubit junction a two-junction 
SQUID, whose critical current can then be tuned via a second flux-bias 
circuit [14,15]. Larger junctions are then permissible, with Ej ~ 10 2 £'c' 
to lO 3 -^- A third method is to fabricate the loop inductance from two or 
more larger critical-current junctions [16]. These junctions are biased with 
phase less than ir/2, and thus act as positive inductors. The advantage to 
this approach is that junction inductors are smaller than physical inductors, 
and fabrication imperfections in the critical currents of the junctions tend 
to cancel out and make the tuning of L with Ljq easier. 

In summary, the major difference between the phase, flux, and charge 
qubits is the shape of their nonlinear potentials, which are respectively cu- 
bic, quartic, and cosine. It is impossible at this time to predict which 
qubit type is best because their limitations are not precisely known, es- 
pecially concerning decoherence mechanisms and their scaling. However, 
some general observations can be made. 

First, the flux qubit has the largest nonlinearity. This implies faster 
logic gates since suppressing transitions from the qubit states and 1 to 
state 2 requires long pulses whose time duration scales as 1/ \ojiq — wail [12]. 
The flux qubit allows operation times less than ~ 1 ns, whereas for the 
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phase qubit 10 ns is more typical. We note, however, that this increase in 
speed may not be usable. Generating precise shaped pulses is much more 
difficult on a 1 ns time scale, and transmitting these short pulses to the 
qubit with high fidelity will be more problematic due to reflections or other 
imperfections in the microwave lines. 

Second, the choice between large and small junctions involve tradeoffs. 
Large junctions (Ej 3> Ec) require precise tuning of parameters (L/Lj for 
the flux qubit) or biases (I/Iq for the phase qubit) to produce the required 
nonlinearity. Small junctions [Ej ~ Ec) do not require such careful tuning, 
but become sensitve to 1// charge fluctuations because Ec has relatively 
larger magnitude. 

Along these lines, the coherence of qubits have been compared consider- 
ing the effect of low- frequency 1// fluctuations of the critical current [17]. 
These calculations include the known scaling of the fluctuations with junc- 
tion size and the sensitivity to parameter fluctuations. It is interesting 
that the calculated coherence times for the flux and phase qubits are simi- 
lar. With parameters choosen to give an oscillation frequency of ~ 1 GHz 
for the flux qubit and ~ 10 GHz for the phase qubit, the number of coherent 
logic-gate operations is even approximately the same. 

4 BCS Theory and the Superconducting State 

A more complete understanding of the Josephson effect will require a deriva- 
tion of Eqs. I2.lal and l2.lbl In order to calculate this microscopically, we 
will first review the BCS theory of superconductivity [18] using a "pair 
spin" derivation that we believe is more physically clear than the standard 
energy-variational method. Although the calculation follows closely that 
of Anderson [19] and Kittel [20], we have expanded it slightly to describe 
the physics of the superconducting phase, as appropriate for understanding 
Josephson qubits. 

In a conventional superconductor, the attractive interaction that pro- 
duces superconductivity comes from the scattering of electrons and phonons. 
As illustrated in Fig. |3fa), to first order the phonon interaction scatters an 
electron from one momentum state to another. When taken to second or- 
der (Fig. Elb)), the scattering of a virtual phonon produces a net attractive 
interaction between two pairs of electrons. The first-order phonon scatter- 
ing rates are generally small, not because of the phonon matrix element, 
but because phase space is small for the final electron state. This implies 
that the energy of the second order interaction can be significant if there 
are large phase-space factors. 

The electron pairs have the largest net interaction if every pair is allowed 
by phase space factors to interact with every other pair. This is explicitly 
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Fig. 3. Feynman diagram of electron-phonon interaction showing (a) first- and 
(b) second-order processes. 

created in the BCS wavefunction by including only pair states (Cooper 
pairs) with zero net momentum. Under this assumption and using a second 
quantized notation where ct is the usual creation operator for an electron 
state of wavevector k, the most general form for the electronic wavefunction 
is 

y = H(u k + v k e^c{cl k )\0) , (4.1) 
k 

where u k and v k are real and correspond respectively to the probability 
amplitude for a pair state to be empty or filled, and are normalized by u\+ 
Vp. = 1. For generality we have included a separate phase factor <fi k for 
each pair. Because each pair state is described as a two state system, the 
wavefunction may also be described equivalently with a "pair-spin" tensor 
product 

-n(,£>. <"> 

and the Hamiltonian given with Pauli matrices a x k, &yk, and <r z k- 

The kinetic part of the Hamiltonian must give \& in the ground state 
with pairs occupied only for \k\ < fc/, where fc/ is the Fermi momentum. 
If we define the kinetic energy of a single electron, relative to the Fermi 
energy, as £ fc , then the kinetic Hamiltonian for the pair state is 

HK = -Y,Sk°zk- (4-3) 

The solution of Hk^ — E k ± i S gives for the lowest energy, Ek—, the values 
Vk = 1 for |fc| < fc/, and v k — for |fc| > fc/, as required. An energy 
Ek+ — Ek- = 2 |£ fe | is needed for the excitation of pairs above the Fermi 
energy or the excitation of holes (removal of pairs) below the Fermi energy. 

The potential part of the pair-spin Hamiltonian comes from the second- 
order phonon interaction that both creates and destroys a pair, as illustrated 
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in Fig. GUb). The Hamiltonian for this interaction is given by 



(4.4) 



k,l 



and can be checked to correspond to the second-quantization Hamiltonian 
Ha = — V c k c -k c k c -k by using the translation a xk — > CfcC_jt + c IcL fc and 



We will first understand the solution to the Hamiltonian Hk + Ha for 
the phase variables <p k . This Hamiltonian describes a bath of spins that are 
all coupled to each other in the x-y plane (Ha) and have a distribution of 
magnetic fields in the z-direction (Hk)- Because Ha is negative, each pair 
of spins becomes aligned with each other in the x-y plane, which implies that 
every spin in the bath has the same phase (f> k . This condition explains why 
the BCS wavefunction has only one phase <j> — (j) k for all Cooper pairs [21]. 
Because there is no preferred direction in the x-y plane, the solution to the 
Hamiltonian is degenerate with respect to <j) and the wavefunction for <f> is 
separable from the rest of the wavefunction. Normally, this means that (j) 
can be treated as a classical variable, as is done for the conventional un- 
derstanding of superconductivity and the Josephson effects. For Josephson 
qubits, where <j> must be treated quantum mechanically, then the behavior 
of 4> is described by an external-circuit Hamiltonian, as was done in Sec. 

For a superconducting circuit, where one electrode is biased with a volt- 
age V, the voltage can be accounted for with a gauge transformation on 
each electron state cl — > e^ 6 /^ ^ vdt c\. The change in the superconducting 
state is thus given by 



(T yk — > i(c k c-k - c fe c -fc)- 



l[(u k + Vke^eWV / Vdt c\ I ™* c t_ fe ) | > (4.5) 



HK + Vke^WK) / ™*] c t c t_ fc) | ) . (4.6) 



fc 



The change in 



can be written equivalently as 



# 2eV 



(4.7) 



dt h 



which leads to the AC Josephson effect. 

The solution for Uk and v k proceeds using the standard method of mean- 
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Fig. 4. Bloch sphere solution of the Hamiltonian (a x ,a y ,a z ) • (B x , B y , B z ). The 
vector B gives the direction of the positive energy eigenstate. 



field theory, with 



V x - 

( H a) = 2_JyOxk Wxl) + <?yk ((Tyl)) , (4.8) 

k,l 

{* xl ) = {u u v l e- i *).a x .( Kv ^ (4.9) 

= 2uivi cos0 , (4-10) 
(<T y i) = 2uivi sin0 . (4-11) 

Using the standard definition of the gap potential, one finds 

A = VJ2uiv l , (4.12a) 
i 

H = H K + (H A ) (4.12b) 

= - y^((T x k,^yk,(Tzk) ■ (A cos 0, A sin 0,4) . (4.12c) 



This Hamiltonian is equivalent to a spin 1/2 particle in a magnetic field, 
and its solution is well known. The energy eigenvalues of Hty = Ek±^ are 
given by the total length of the field vector, 

£ fc± =±(A 2 +^) 1/2 , (4.13) 



and the directions of the Bloch vectors describing the Ek+ and Ek- eigen- 
states are respectively parallel and antiparallel to the direction of the field 
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vector, as illustrated in Fig. 0] The ground state solution ^>k- is given by 

u k = Jl(l + ^-), (4.14) 




'•a =w'-( 1 -tH ' ( 4 - 15 ) 

(4.16) 

with the last equation required for consistency. The excited state is 
similarly described, but with Uk and Vk interchanged and <f> — > + 7r. 

At temperature T = the energy gap A may be solved by inserting the 
solutions for Uk and Vk into Eq. I4.12al 

A = yV A „ - . (4.17) 

V 2 ( A2 +^) 1/2 

Converting to an integral by defining a density of states Nq at the Fermi 
energy, and introducing a cutoff of the interaction V at the Debye energy 
8d, one finds the standard BCS result, 

A = 26 D e- 1/NoV . (4.18) 

Two eigenstates E^- and E^ + have been determined for the pair Hamil- 
tonian. Two additional "quasiparticle" eigenstates must exist, which clearly 
have to be single-particle states. These states may be solved for using di- 
agonalization techniques, giving 

*fcO = 4|0) , (4.19a) 
*fci = cL fc |0) . (4.19b) 



Fortunately, these states may be easily checked by inspection. The kinetic 
part of the Hamiltonian gives Hn^m,! = since &kO,i corresponds to 
the creation of an electron and a hole, and the electron-pair and hole-pair 
states have opposite kinetic energy. The potential part of the energy also 
gives (-Ha) ^fco.i = since the interaction Hamiltonian scatters pair states. 
Thus the eigenenergies of *feo,i are zero, and these states have an energy 
Ek = |-Efe-| above the ground state. 

The quasiparticle operators that take the ground-state wavefunction to 
the excited states are 

7fc = u k c l - Vke-^c-k , (4.20) 
ill = u fe c -fe + Vke-'+Ck, (4.21) 
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Fig. 5. Energy- level diagram for the ground-pair state (solid line), two quasipar- 
ticle states (dashed lines), and the excited-pair state (short dashed line). 

which can be easily checked to give 



A summary of these results is illustrated in Fig. where we show 
the energy levels, wavef unctions, and operators for transitions between the 
four states. The quasiparticle raising and lowering operators 7^, 7L, 
7 fe0 , and j kl produce transitions between the states and have orthogonality 
relationships similar to those of the electron operators. 

It is interesting to note that the ground and excited pair states are con- 
nected by the two quasiparticle operators e~ % ^ikiiko^ k- = ^k+- Because 
the value u;i>/ changes sign between ^k- and ^k+- an d is zero for ^Sko,ii 
the gap eauation l4.12al including the effect of quasiparticles is proportional 
to (1 — 7^, 7fco — 7fei7fci)- Along with the energy levels, these results imply 
that the two types of quasiparticles are independent excitations. 

5 The Josephson Effect, Derived from Perturbation Theory 

We will now calculate the quasiparticle and Josephson current for a tun- 
nel junction using hrst and second order perturbation theory, respectively. 
We note that our prior calculations have not been concerned with electrical 
transport. In fact, the electron operators describing the superconduct- 
ing state have not been influenced by charge, and thus they correspond to 
the occupation of an effectively neutral state. Because a tunneling event 
involves a real transfer of an electron, charge must now be accounted for 



7L(«* + «fc^4cL fc )|0) 

7k(«*+«* e< *4 c U)l°) 




(4.22) 
(4.23) 
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properly. We will continue to use electron operators for describing the 
states, but will keep track of the charge transfer separately. 

When an electron tunnels through the barrier, an electron and hole state 
is created on the opposite (left and right) side of the barrier. The tunneling 
Hamiltonian for this process can be written as 

H T = H T + + H T - + H T+ + H T - (5.1a) 

= 0'LR C L C R + t_ L _ R C_ L C [ _ R + t* LR C^ L C R + t*_ L _ R C^_ L C_ R ^j.l})) 

L.R 

where t LR is the tunneling matrix element, and the L and R indices refer 

respectively to momentum states k on the left and right superconductor. 
— > — > 

The first two terms Ht+ and Ht- correspond to the tunneling of one 
electron from left to the right, whereas Ht+ and Ht- are for tunneling 
to the left. The Hamiltonian is explicitly broken up into Ht+ and Ht- 
to account for the different electron operators c\. and cL fc for positive and 
negative momentum. 

The electron operators must first be expressed in terms of the quasipar- 
ticle operators 7 because these produce transitions between eigenstates of 
the superconducting Hamiltonian. Equations 14. 201 14.211 and their adjoints 
are used to solve for the four electron operators 

Cfe = Ufc7feo + u fe e# 7fci c-k = Wfe7fei ~ Vfce # 7L r 5 2 ) 



Substituting Eqs. 15.21 into I5.1bl one sees that all four terms of the 
Hamiltonian have operators 7^ that produce quasiparticles. We calculate 

— y 

here to first order the quasiparticle current from L to R given by Ht+ + 
Ht-- The Feynman diagrams (a) and (b) in Fig. respectively describe 
the tunneling Hamiltonian for the Ht+ and Ht- terms. In this diagram 
a solid line represents a Cooper pair state in the ground state, whereas a 
quasiparticle state is given by a dashed line. Only one pair participates in 
the tunneling interaction, so only one of the three solid lines is converted 
to a dashed line. The line of triangles represents the tunneling event and 
is labeled with its corresponding Ht Hamiltonian, with the direction of 
the triangles indicating the direction of the electron tunneling. The c k 
operators, acting on the L or R lead, is rewritten in terms of the 7 operators 
and placed above or below the vertices. Since only y< operators give a 
nonzero term when acting on the ground state, the effect of the interaction is 
to produce final states ^/j' R with total energy Er+El, and with amplitudes 
given at the right of the figure. 

The two final states in FigUJa) and (b) are orthogonal, as well as states 
involving different values of L and R. The total current is calculated as 
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Fig. 6. First-order Feynmann diagrams for interaction Ht+ (a) and Ht- (b). 
Solid lines are Cooper-pair states, dashed lines are quasiparticle excitations, and 
arrow-lines represents tunneling interaction. Electron operators arising from 
interaction are displayed next to vertices. 

an incoherent sum over all possible final quasiparticle states, under the 
condition that the total quasiparticle energy for the final state is equal to 
the energy gained by the tunneling of the electron 



The total current from L to R is given by e multiplied by the transition rate 



E R + E L = eV . 



(5.3) 




(E R +E L =eV) 






27re 



L,R 
(E R +E L =eV) 





-L-R 



{ U R V L? 



(5.4b) 



L,R 




ErX^Ac) 
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where in the last equation we have expressed the conservation of energy 
with a Dirac <5-function, and have assumed matrix elements \t\ 2 of constant 
strength. Because E(£ k ) = E(—£ k ) and u k (£,k) = v k (— £ fe ), one finds 

oc oo 

V = |t| 2 iVoflJVoL y («i + ui)d£ L J (u 2 R + v 2 R ) d£ R S(eV - E L - E R ) 



oo oo 

= \t\ 2 N 0R N 0L J di L J dt R S(eV - E L - E R ) . 



(5.5a) 



o o 



This result is equivalent to the standard "semiconductor model" of the 
quasiparticle current, which predicts no current for V < 2A/e, a rapid rise 
of current at 2A/e, and then a current proportional to V at large voltages. 
Note that Eq. I5.4cl has a sum over the occupation probability v\ of the pair 
state and the occupation probability u R of a hole-pair state, as is expected 
given the operators c L c R in the tunneling Hamiltonian. The final result of 
Eq. !5.5al does not have these factors because the occupation probability is 
unity when summed over the ±£ fc states. 

It is convenient to express the tunneling matrix element in terms of the 
normal-state resistance of the junction, obtained by setting A = 0, with the 
equation 

l/R N ee T2/V (5.6a) 



47re 
Aire 2 



LXJ oo 

\t\ 2 N 0R N QL J d£ L J dt R 6(eV-Z L -£ R )/V(5.&b) 



t\ 2 N 0R N 0L . (5.6c) 



We now calculate the tunneling current with second-order perturbation 
theory. The tunneling Hamiltonian, taken to second order, gives 

H {2) =J2H T -H T , (5.7) 

i £i 

where is the energy of the intermediate state i. Because the terms in 
have both and 7 operators, the second-order Hamiltonian gives a nonzero 
expectation value for the ground state. This is unlike the first-order theory, 
which produces current only through the real creation of quasiparticles. 

Because Ht has terms that transfer charge in both directions, HtHt 
will produce terms which transfer two electrons to the right, two to the left, 
and with no net transfer. With no transfer, a calculation of the second- 
order energy gives a constant value, which has no physical effect. We first 
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Fig. 7. Second-order Feynman diagrams for the transfer of two electrons across 
the junction. Only nonzero operators are displayed next to vertices. 



calculate terms for transfer to the right from (Ht+ + Ht-){Ht+ + Ht-), 

— > — > — > — ► 

which gives nonzero expectation values only for Ht+Ht- + Ht-Ht+- 
The Feynman diagrams for these two terms are given in Fig. Ufa) and (b), 
where we have displayed only the amplitudes from the nonzero operators. 
The expectation value of these two Hamiltonian terms is given by 



4 2) = -£<*-K*-lfe 



LM 



LR 



Er + El 

x{u R v L e i<t>L ) I* 1 :) \^ R ) (5.8a) 

-2 If ~ (5.8b) 
1% Er + El 

OC oc 

-2\tfe-N 0R N 0L J d U J ^A^-^- (5.8c) 



H ^ e iS I d6 R [ d9 L — r^r- (5.8d) 

2ire 2 R N J J cosh 6 R + cosh 6 L V ; 



— OO — OO 
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where we have used t* LR = t_ L _ R and assumed the same gap A for both 
superconductors. A similar calculation for transfer to the left gives the 
complex conjugate of Eq. I5.8el The sum of these two energies gives the 
Josephson energy Uj, and using Eq. 12.51 the Josephson current Ij, 

h = |-|- S inS, (5.10) 

where Rk = h/e 2 is the resistance quantum. Equation 15. 101 is the stan- 
dard Ambegaokar-Baratoff formula [22] for the Josephson current at zero 
temperature. 

The Josephson current is a dissipationless current because it arises from 
a new ground state of the two superconductors produced by the tunneling 
interaction. This behavior is in contrast with quasiparticle tunneling, which 
is dissipative because it produces excitations. It is perhaps surprising that a 
new ground state can produce charge transfer through the junction. This is 
possible only because the virtual quasiparticle excitations are both electrons 
and holes: the electron-part tunnels first through the junction, then the 
hole-part tunnels back. Only states of energy A around the Fermi energy 
are both electron- and hole- like, as weighted by the (v > R u R )(u L v L ) term in 
the integral. 

The form of the Josephson Hamiltonian can be understood readily by 
noting that the second-order Hamiltonian, 

Ht+Ht- ~ \t\ 2 Y,c L c_ L c R cl R (5.11) 

L,R 

\t\ 2 
1 L,R 

corresponds to the pair-scattering Hamiltonian of Eq. 14.41 Comparing 
with the gap-equation solution, one expects Uj ~ |i| 2 Acos(5, where the 
cos 5 term arises from the spin-spin interaction in the x-y plane. 

We would like to make a final comment on a similarity between the 
BCS theory and the Josephson effect. In both of these derivations we 
see that a dissipative process that is described in first-order perturbation 
theory, such as phonon scattering or quasiparticle tunneling, produces in 
second order a new collective superfluid behavior. This collective behavior 
emerges from a virtual excitation of the dissipative process. Dissipation 
is normally considered undesirable, but by designing systems to maximize 
dissipation, it may be possible to discover new quantum collective behavior. 

With this understanding of the Josephson effect and quasiparticle tun- 
neling, how accurate is the description of the Josephson junction with the 
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Hamiltonian corresponding to Eq. 15. 9V There are several issues that need 
to be considered. 

First, quasiparticle tunneling is a dissipative mechanism that produces 
decoherence. Although it is predicted to be absent for V < 2A/e, mea- 
surements of real junctions show a small subgap current. This current is 
understood to arise from multiple Andreev reflections, which are described 
as higher-order tunneling processes. We thus need a description of the 
tunnel junction that easily predicts these processes for arbitrary tunneling 
matrix elements. This is especially needed as real tunnel junctions do 
not have constant matrix elements, as assumed above. Additionally, we 
would like to know whether a small number of major imperfections, such as 
"pinhole" defects, will strongly degrade the coherence of the qubit. 

Second, quasiparticle tunneling has been predicted for an arbitrary DC 
voltage across the junction. However, the qubit state has (V) = 0, but 
may excite quasiparticles with AC voltage fluctuations. This situation is 
difficult to calculate with perturbation theory. In addition, is it valid to 
estimate decoherence from quasiparticles at zero voltage simply from the 
junction resistance at subgap voltages? 

Third, how will the Josephson effect and the qubit Hamiltonian be mod- 
ified under this more realistic description of the tunnel junction? 

All of these questions and difficulties arise because perturbation theory 
has been used to describe the ground state of the Josephson junction. The 
BCS theory gives basis states that best describe quasiparticle tunneling for 
large voltages, not for V — > 0. A theory is needed that solves for the 
Josephson effect exactly, with this solution then providing the basis states 
for understanding quasiparticle tunneling around V = 0. This goal is 
fulfilled by the theory of quasiparticle bound states, which we will describe 
next. 

6 The Josephson Effect, Derived from Quasiparticle Bound States 

We begin our derivation of an exact solution for the Josephson effect with an 
extremely powerful idea from mesoscopic physics: electrical transport can 
be calculated under very general conditions by summing the current from a 
number of independent "conduction channels" , with the transport physics 
of each conduction channel determined only by its channel transmission 
probability r< [23,24]. For a Josephson junction, the total junction current 
Ij can be written as a sum over all channels i 

h = Y,h(Ti), (6.i) 

i 

where Ij(r) is the current for a single channel of transmission r, which may 
be solved for theoretically. For a tunnel junction, the number of chan- 
nels is estimated as the junction area divided by the channel area (A//2) 2 , 
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Fig. 8. Plot of potential vs. coordinate x with a positive delta-function tunnel 
barrier VqS(x). Scattering of plane wave states is shown in (a), whereas (b) is a 
plot of the bound-state wavefunction. The delta-function barrier is negative in 
(b), as required for producing a bound state. 

where A/ is the Fermi wavelength of the electrons. Of course, the difficulty 
of determining the distribution of the channel transmissions still remains. 
This often may be estimated from transport properties, and under some 
situations can be predicted from theory [25-27]. 

Because transport physics is determined only by scattering parameter- 
ized by t, we may make two simplifying assumptions: the transport can be 
solved for using plane waves, and the scattering from the tunnel junction 
can be described by a delta function. The general theory has thus been 
transformed into the problem of one-dimensional scattering from a delta 
function, and an exact solution can be found by using a simple and clear 
physical picture. 

Central to understanding the Josephson effect will be the quasiparticle 
bound state. To understand how to calculate a bound state [28,29], we will 
first consider a normal-metal tunnel junction and with a (^-function barrier 
Vq5(x), as illustrated in Fig. |H1 For an electron of mass m and wavevector 
k, the wavefunctions on the left and right side of the barrier are 



where A, B, and C are respectively the incident, reflected, and transmitted 
electron amplitudes. From the continuity equations 



* L = Ae lkx + Be 
9 R = Ce ikx , 



(6.2) 
(6.3) 



*L = <Fr 

d^ L d^ R 2mV 



(6.4) 



dx dx h 2 

evaluated at x = 0, the amplitudes are related by 



(6.5) 



A + B 



C 



(6.6) 



ikA - ikB = ikC - 



2mV 



c 



(6.7) 
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Fig. 9. Plot of quasiparticle energies E K verses momentum k near the ±ft/ Fermi 
surfaces. The two-component eigenfunctions are also displayed for each of the 
four energy bands. Also indicated are the quasiparticle states A-E used for the 
bound-state calculation. 



The transmission amplitude and the probability are 

C 1 



A 1 + irj 

C 2 
A 



(6.8) 



i + n 



Y — , (6.9) 



where r\ — mVo/h 2 k. The bound state can be determined by finding the 
pole in the transmission amplitude. A pole describes how a state of finite 
amplitude may be formed around the scattering site with zero amplitude of 
the incident wavefunction, which is the definition of a bound state. The pole 
at 77 = i gives kb — —imVo/h 2 , and a wavefunction around the scattering 
site = Ce( mV °/ h This describes a bound state only when Vo is 
negative, as expected. 

A superconducting tunnel junction will also have bound states of quasi- 
particles excitations. These bound states describe the Josephson effect 
since virtual quasiparticle tunneling was necessary for the perturbation cal- 
culation in the last section. The Bogoliubov-deGcnnes equations describe 
the spatial wavefunctions, whose eigenstates are given by the solution of the 
Hamiltonian 

Htp ± e iKX = {±k f ^a z + A<7 x ^j ^e lKX , (6.10) 

where ( / 9 ± e 4KX are the slowly varying spatial amplitudes of the exact wave- 
function ip ± e lKX e ±lkfX . As illustrated in Fig. the ±fc/K term corresponds 
to the kinetic energy at the ±fc/ Fermi surfaces using the approximation 



Martinis et al: Superconducting Qubits ... 23 



(kf + «) 2 /2 ~ const. + kfK. As expected for a spin-type Hamiltonian, the 
two eigenvalues are 

£ K = ±(^ + A 2 ) 1/2 , (6.H) 

where £ K = H 2 kfK/m is the kinetic energy of the quasiparticle referred to 
the Fermi energy. The eigenvectors are also displayed in Fig. [§1 where u K 
and v K are given by Eqs. 14.141 and 14.151 Because the two energy bands 
represent quasiparticle excitations, the lower band is normally filled and its 
excitations correspond to the creation of hole states. 

We can solve for the quasiparticle bound states by first writing down the 
scattering wavefunctions in the left and right superconducting electrodes. 
An incoming quasiparticle state, point A in Fig. is reflected off the tunnel 
barrier to states B and C and is transmitted to states D and E [30]. The 
wavefunctions are then given by 

* L = a( U .. )e iKX +B[ V . )e iKX + c( " \e~ iKX (6.12) 

- d( V j \c- ikx + e( U , }e iKX , (6.13) 

where we have used the relations v = v K = u_ K and u = u K = w_ K , and 
we have included the phases <f> L and <f) R of the two states. The continuity 
conditions Eqs. 16.41 and 16.51 solved for both the components of the spin 
wavefunction, gives the matrix equation 
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ue 
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(6.14) 



The scattering amplitudes for B-E have poles given by the solution of 



{u 4 + v 4 )(l + rf) - 2{uvf{rf + cos<5) 

,2 i „,2 







Using the relations u 2 + v 2 — 1, Ej = E^ = A/2uv, and r 
the energies of the quasiparticle bound states are 



E J± =±A[l-rsin 2 (5/2)] 



1/2 



(6.15) 
1/(1 + v 2 ), 

(6.16) 



Because these two states have energies less than the gap energy A, they are 
energetically "bound" to the junction and thus have wavefunctions that are 
localized around the junction. 

The dependence of the quasiparticle bound-state energies on junction 
phase is plotted in Fig. ^]for several values of r. The ground state is nor- 
mally filled, similar to the filling of quasiparticle states of negative energy. 
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Fig. 10. Plot of quasiparticle bound-state energies Ej- and Ej+ vs. the phase dif- 
ference 8 across the junction, for three values of tunneling transmission r. Quasi- 
particles are produced by vertical transitions from the E.j- to Ej+ band. As 
indicated by the arrow, the energy gap Ej+ — Ej- is always greater than V2A 
at 8 = tt/2. 



The energy Ej- corresponds to the Josephson energy, as can be checked in 
the limit r — > to give 

At At 

Ej^~-A + — —cosS. (6.17) 

This result is equivalent to Eq. 15.91 after noting that the normal-state con- 
ductance of a single channel is 1/Rn — 2t/Rk- 

The current of each bound state is given by the derivative of its energy 

j ^dE J± 

Ij± = w ^r^ (6 - 18) 

in accord with Eq. 12.51 Since the curvature of the upper band is opposite 
to that of the lower band, the currents of the two bands have opposite 
sign Ij + = —Ij-. For level populations of the two states given by f±, 
the average Josephson current is (Ij) = //_(/_ — / + ). For a thermal 
population, f± are given by Fermi distributions, and the Josephson current 
in the tunnel junction limit gives the expected Ambegaokar-Baratoff result 

Uj) = l^p. s i n s.( — L —1 ^) (6.19) 

w/ 2 eR N \e- A / kT + l e A / kT + lJ V ; 

_ 7rA(r) tanh ( A / 2 fcr)sin5 . (6.20) 
2 ei?jv 
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Fig. 11. Plot of semiclassical solutions for the tunneling through a barrier (a) 
and tunneling through an energy gap (b). Imaginary solutions to k and 8 are 
used to calculate the tunneling rates. 

7 Generation of Quasiparticles from Nonadiabatic Transitions 

In this description of the Josephson junction, the Josephson effect arises 
from a quasiparticle bound state at the junction. Two bound states ex- 
ist and have energies Ej + and Ej~, with the Josephson current from the 
excited state being of opposite sign from that of the ground state. We 
will discuss here the small- voltage limit [31,32], which can be fully under- 
stood within a semiclassical picture by considering that a linear increase in 
5 produces nonadiabatic transitions between the two states. 

The junction creates "free quasiparticles" , those with E > A, via a two- 
step process. First, a transition is made from the ground to the excited 
bound state. This typically occurs because a voltage is placed across the 
junction, and the linear change of 5 causes the ground state not to adiabat- 
ically stay in that state. For a high-transmission channel, the transition is 
usually made around 5 « 7r, where the energy difference between the states 
is the lowest and the band bending is the highest. Because this excited 
state initially has energy less than A, the state remains bound until the 
phase changes to 27r and the energy of the quasiparticle is large enough to 
become unbound and diffuse away from the junction. The quasiparticle 
generation rate is thus governed by dS/dt and will increase as V increases. 

The quasiparticle transition rate can be predicted using a simple semi- 
classical method. We will first review WKB tunneling in order to later 
generalize this calculation to energy tunnelling. In Fig. Illf a). we plot a 
cubic potential V(x) versus x and its solution k 2 = 2m[E — V(x)]/h 2 . The 
solution for k is real or imaginary depending on whether E is greater or 
less than V(x). A semi-classical description of the system is the particle 
oscillating in the well, as described by the loop in the solution of Re{fc}. A 
solution in the imaginary part of k connects a turning point on this loop, 
labeled A, with the turning point of the free-running solution, labeled B. 
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The probability of tunneling each time the trajectory passes point A is given 
by the standard WKB integral of the imaginary action 



W 

s 



= exp[-25] 



(l/h) 



dx Imp 



dxlm k 



(7.1) 
(7.2) 

(7.3) 



The transition rate for a nonadiabatic change in a state may be cal- 
culated in a similar fashion. In Fig. II If we plot the solution of Eq. 
16.161 for 8 versus E. In the "forbidden" region of energy \E\ < A-y/l — r, 
the solution of 5 has an imaginary component. As the bias of the system 
changes and the system trajectory moves past point A, then this state can 
tunnel to point B via the connecting path in the imaginary part of 6. The 
probability for this event is given by Eq. 17.11 with S given by the integral 
of the imaginary action 



S=(l/h) J dEt imag 



where we define an imaginary time by 



Im 5 



dS/dt 



(7.4) 



(7.5) 



Rewriting Eq. 16.161 as (E/A) 2 = 1 — r(l — cos(5)/2 and using dS/dt = 
(2e/H)V, one finds the action is given by the integral 



5 



2eV 
A 

eV ' 



de Im {- arccos [l + (e 2 - l)2/r] } (7.6) 

(1-t)tt/2 (t-^1) 
ln(2/r) (r ->■()) 

(l-r)[ln(2/T) + v^(7r/2-ln2)] , (interp.) 



where the last interpolation formula approximates well a numerical solution 
of Eq. 17.61 The limiting expression for r — ► 1 gives the standard Landau- 
Zener formula appropriate for a two-state system. In the tunnel-junction 
limit t — > one finds 

W = (r/2) 2A / eV . (7.7) 

For the case of a constant DC bias voltage V, the total junction current 
(Ij) may be calculated with this transition rate and an attempt rate T = 
(2e/h)V given by the frequency at which 5 passes n/2. Using Eq. 17.61 and 
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Fig. 12. (a) Plot of average junction current (Ij) versus inverse DC voltage V for 
transmission coefficients r = 0.8, 0.5, 0.2, and O.Of (from Ref. [32] ). Solid lines 
are from exact calculation, and dashed lines are from predictons of Eqs. 17.81 and 
17.61 The time dependence of the Josephson current Ij is plotted for the ground 
state (b) and for a transition (c), where the insets show the trajectory of the 
bound states as Ej vs. S. 



setting the power of quasiparticle generation 2ArW ,r to the electrical power 
(Ij) V, one finds 

(« = ?§W . (7.8) 



This prediction is plotted in Fig. IT2T a1 and shows very good agreement 
with the results of exact calculations [32]. Only the steps in voltage are not 
reproduced, which are understood as arising from the quantization of energy 
eV from multiple Andreev reflection of the quasiparticles. The steps are 
not expected to be reproduced by the semiclassical theory since this theory 
is an expansion around small voltages, or equivalently, large quantization 
numbers. 

The junction current may also be determined from the energies of the 
two bound states. For a constant voltage across the junction, we use Eq. 
12.51 to calculate the charge transferred across the junction after a phase 
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change of 2-k 



Qj = I 

Jo 



2ir/(dS/dt) 



Ijdt 



(7.9) 



2tt r^/Wdt) dUj 



(7.10) 



$o Jo dS 
[Uj(2ir) - Uj(0)} 



(7.11) 



V 



which gives the expected result that the change of energy equals QjV. 
When the junction remains in the ground state, the energy is constant 
Uj(2tt) — Uj(0) = and no net charge flows through the junction. Net 
charge is transferred, however, after a transition. The charge transfer 
2A/V multiplied by the transition rate gives an average current QjTW 
that is equivalent to Eq. 17.81 

Equation 16.181 may be used to calculate the time dependence of the 
Josephson current, as illustrated in Fig. 112( b) and (c). When the sys- 
tem remains in the ground state (b), the junction current is sinusoidal and 
averages to zero. For the case of a transition (c), the current before the 
transition is the same, but the Josephson current remains positive after the 
transition (see Eq. 5.4 of Ref. [32]). The transition itself also produces 
charge transfer from multiple- Andreev reflections (MAR) [31,33] 



This result is perhaps surprising - the junction current at finite voltage arises 
from transfer of charge Qmar cmd a change in the Josephson current. The 
relative contribution of these two currents is determined by the relative size 
of the gap in the bound states. For r — > 1 , all of the junction current is 
produced by Josephson current, whereas for r — > (tunnel junctions) the 
current comes from Qmar- 

For small voltages, the transition event must transfer a large amount of 
charge Qmar in order to overcome the energy gap. In comparing this semi- 
classical theory with the exact MAR theory, Qmar/c has an integer value 
and represents the order of the MAR process and the number of electrons 
that are transferred in the transition. This description is consistent with 
Eq. I7.7l describing the transition probability for an n-th order MAR process, 
where n = 2A/eV, and r/2 represents the matrix element for each order. 

From this example it is clearly incorrect to picture the quasiparticle and 
Josephson current as separate entities, as suggested by the calculations of 
perturbation theory. To do so ignores the fact that quasiparticle tunnel- 
ing, arising from a transition between the bound states, also changes the 
Josephson contribution to the current from 5 = it to 2-7T. 



Qmar = 2A(1-t) 1 / 2 /^ ■ 



(7.12) 
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8 Quasiparticle Bound States and Qubit Coherence 

The quasiparticle bound-state theory can be used to predict both the Joseph-^ 

son and quasiparticle current in the zero-voltage state, as appropriate for 
qubits. In this theory an excitation from the Ej- bound state to the Ej+ 
state is clearly deleterious as it will change the Josephson current, fluctu- 
ating the qubit frequency and producing dccohercnce in the phase of the 
qubit state. For an excitation in one channel, the fractional change in the 
Josephson current is ~ l/7V c h, where iV c h is the number of conduction chan- 
nels. The subgap current-voltage characteristics can be used to estimate 
N c h, which gives an areal density of ~ 10 4 //im 2 [10,34]. For a charge qubit 
with junction area 10~ 2 /im 2 , the qubit frequency changes fractionally by 
~ l/N ch ~ 10~ 2 for a single excitation, and gives strong decoherence. Al- 
though the phase qubit has a smaller change (l/iV c h)Io/4(/o — I) ~ 2x 10~ 5 , 
the excitation of even a single bound state is clearly unwanted. 

Fortunately, these quasiparticle bound states should not be excited in 
tunnel junctions by the dynamical behavior of the qubit. The Sj_ to Ej + 
transition is energetically forbidden because the energy of the qubit states 
are typically choosen to be much less than 2A. Thus, the energy gap of 
the superconductor protects the qubit from quasiparticle decoherence. 

If a junction has "pinhole" defects, where a few channels have r — > 1, 
then the energy gap will shrink to zero at S = 7T. However, only the flux 
qubit will be sensitive to quasiparticles produced at these defects since it 
operates near 5 = ir. In contrast, the phase qubit always retains an energy 
gap of at least y/2A around its operating point 6 — tt/2 (see arrow in Fig. 
110(1 . We note this idea implies that a phase qubit can even be constructed 
from a microbridge junction, which has some channels [26] with r = 1. 
Although the phase qubit is completely insensitive to pinhole defects, this 
advantage is probably unimportant because Al-based tunnel junctions have 
oxide barriers of good quality. 

Pinhole defects also change the Josephson potential away from the -cos S 
form. This modification is typically unimportant because the deviation is 
smooth and can be accounted for by a small effective change in the critical 
current. 

The concept that the energy gap A protects the junction from quasi- 
particle transitions suggests that superconductors with nonuniform gaps 
may not be suitable for qubits. Besides the obvious problem of conduc- 
tion channels with zero gap, channels with a reduced gap may cause stray 
quasiparticles to be trapped at the junction. The high-T c superconductors, 
with the gap suppressed to zero at certain crystal angles, are an obvious 
undesirable candidate. However, even Nb could be problematic since it has 
several oxides that have reduced or even zero gap. Nb based tri-layers may 
also be undesirable since the thin Al layer near the junction slightly reduces 
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the gap around the junction. In contrast, Al may not have this difficulty 
since its gap increases with the incorporation of oxygen or other scattering 
defects. It is possible that these ideas explain why Nb-based qubits do not 
have coherence times as long as Al qubits [6, 10]. 



9 Summary 

In summary, Josephson qubits are nonlinear resonators whose critical el- 
ement is the nonlinear inductance of the Josephson junction. The three 
types of superconducting qubits, phase, flux, and charge, use this nonlin- 
earity differently and produce qubit states from a cubic, quartic, and cosine 
potential, respectively. 

To understand the origin and properties of the Josephson effect, we have 
first reviewed the BCS theory of superconductivity. The superconducting 
phase was explicitly shown to be a macroscopic property of the supercon- 
ductor, whose classical and quantum behavior is determined by the external 
electrical circuit. After a review of quasiparticle and Josephson tunneling, 
we argued that a proper microscopic understanding of the junction could 
arise only from an exact solution of the Josephson effect. 

This exact solution was derived by use of mesoscopic theory and quasi- 
particle bound states, where we showed that Josephson and quasiparticle 
tunneling can be understood from the energy of the bound states and their 
transitions, respectively. A semiclassical theory was used to calculate the 
transition rate for a finite DC voltage, with the predictions matching well 
that obtained from exact methods. 

This picture of the Josephson junction allows a proper understanding 
of the Josephson qubit state. We argue that the gap of the supercon- 
ductor strongly protects the junction from quasiparticle tunneling and its 
decoherence. We caution that an improper choice of materials might give 
decoherence from quasiparticles that are trapped at sites near the junction. 

We believe a key to future success is understanding and improving this 
remarkable nonlinearity of the Josephson inductance. We hope that the 
picture given here of the Josephson effect will help researchers in their quest 
to make better superconducting qubits. 
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